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S ■ Abstract 
S3 . 

Q , We propose a theory of two-magnon Raman scattering from the insulating 

O 

parent compounds of high-T c superconductors, which contains information 
not only on magnetism, but also on the electronic properties in these materi- 
als. We use spin density wave formalism for the Hubbard model, and study 
diagrammatically the profile of the two-magnon scattering and its intensity 
dependence on the incoming photon frequency both for Uj <C U and in the 
resonant regime, in which the energy of the incident photon is close to the gap 
between conduction and valence bands. In the nonresonant case, we identify 
the diagrams which contribute to the conventional Loudon-Fleury Hamilto- 
nian. In the resonant regime, where most of the experiments have been done, 
we find that the dominant contribution to Raman intensity comes from a dif- 
ferent diagram, one which allows for a simultaneous vanishing of all three of 



its denominators (i.e., a triple resonance). We study this diagram in detail 
and show that the triple resonance, combined with the spin-density-wave dis- 
persion relation for the carriers, explains the unusual features found in the 
two-magnon profile and in the two-magnon peak intensity dependence on the 
incoming photon frequency. In particular, our theory predicts a maximum of 
the two-magnon peak intensity right at the upper edge of the features in the 
optical data, which has been one of the key experimental puzzles. 
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I. INTRODUCTION 



There is a widespread belief that strong electron-electron correlations in the high-T c 
compounds may hold a clue to the phenomenon of high-temperature superconductivity.0~i 
One of the manifestations of these correlations is in the fact that the insulating parent com- 
pounds are antiferromagnets. An important probe of antiferromagnetism is magnetic Raman 
scattering.0~0 Its prominent signature in the underdoped high-T c materials is a two-magnon 
peak observed at about 3000 cmT 1 . To first approximation, this peak can be attributed to 
inelastic scattering from the two-magnon excitations.@ In fact, Raman experiments yielded 
the first estimate of the exchange interaction constant in La^CuO^. The two-magnon 
peak has also been observed in the electron-doped materialsl! and more recently^ in the 
underdoped YBa 2 Cu 3 O e+x materials (up to x = 0.9). 

The traditional framework for understanding the two-magnon Raman scattering in anti- 
ferromagnets has been an effective Hamiltonian for the interaction of light with spin degrees 
of freedom known as the Loudon-Fleury Hamiltonian0 

H=a J E(e i .R ij )(e r R ij )S i .S j . (1) 

<ij> 

Here e; and e/ are the polarization vectors of the in- and outgoing photons, a is the (generally 
poorly known) coupling constant, and Rjj is a vector along the bond connecting two nearest 
neighbor sites i and j. This equation has been reexamined in recent years to better account 
for a possible importance of quantum fluctuations, the four-magnon processes, and the 
further neighbor terms 

The Loudon-Fleury theory is expected to work well when the frequencies of the incoming 
and outgoing photons are considerably smaller than the gap between conduction and valence 
bands. The experimental reality in high-T c materials is such, however, that the two-magnon 
scattering is observed mostly in or near the so-called resonant regime, when the frequencies of 
the ingoing and/or outgoing photons are close to the gap value. This is simply a consequence 
of the fact that the experiments are mostly done with the visible or near ultraviolet lasers, 
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while the gap and the upper Hubbard band states lie in the same range of frequences 
(roughly between 2 eV and 3 eV) in the high-T c materials. The experimental cross-sections 
vary strongly in this range of incident photon frequencies,0'0'0 and it becomes an issue 
whether the Loudon-Fleury Hamiltonian is still applicable. Whatever the answer to this 
question, we still need to have the means of evaluating the variation of the overall scale of 
the coupling constant a in Eq. (HD with the incident photon frequency uj{. 

The profile of the Raman cross-section as a function of the transferred photon frequency 
is shown in Fig. [l|. The two-magnon peak is clearly identifiable. The behavior of the two- 
magnon peak height as a function of the incident photon frequency is presented in Fig. |2|, 
where the absorptive part of the dielectric constant is also shown for comparison. 

The key experimental features that require explanation are the following: 

In Fig. |: 

(a) Asymmetry of the two-magnon peak profile: the two-magnon peak is asymmetric, with 
the spectral weight shifted to higher frequencies. 

(b) Selection rules: Loudon-Fleury Hamiltonian predicts no scattering in the A\ g configu- 
ration; experimentally in the resonant regime the A\ g cross-section is about half of that in 
the Bi g geometry. 

(c) Stability of the two-magnon peak profile: as the incident photon frequency changes over 
the frequency range from 15, 000 to 22, 000 cm -1 , the two-magnon scattering intensity profile 
merely scales with incident frequency, without a noticeable change in shape. (However, at 
higher photon energies, a distortion of the Raman spectrum occurs.) 

In Fig. |: 

(a) A single peak: ordinarily one might expect two peaks, the so called ingoing and outgoing 
resonances, EIE and here only one is observed. 

(b) Peak location: a comparison with the dielectric constant shows that the strength of the 
two-magnon Raman scattering is at its maximum away from the band edge, in fact right at 
the upper end of the features in the optical data that can be interpreted as the particle-hole 
excitations between the lower and upper Hubbard bands. (Phonons, on the contrary, are 
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known to resonate with the features in the dielectric constant. 

Raman scattering in a one-band model of correlated electrons has been considered several 
times in the literature .Bill! Shastry and Shraimanil have recently given a derivation of the 
Loudon-Fleury Hamiltonian starting from the usual large-?/ Hubbard model. Working in a 
localized basis, they performed a hopping expansion controlled by t/(U — to), where t and U 
are the nearest-neighbor hopping and on-site Coulomb repulsion, and u is of the order of the 
photon frequencies. The leading term in the expansion turned out to be the Loudon-Fleury 
Hamiltonian, Eq. ([]]). 

A simple, though somewhat crude interpretation of this outcome is to envision spins in 
a classical Neel state. The incoming photon is absorbed by moving, say, a down spin to the 
neighboring site occupied by an up spin. The up spin then returns to the formerly down 
site, emitting the outgoing photon. The overall result is light scattering accompanied by 
the nearest neighbor spin exchange, which is encoded in Eq. (|l|). This process is shown in 
Fig. |. 

Notice that the intermediate state in this process has a doubly occupied site, and is thus 
of order U above the ground state in energy. When the photon energy is far smaller than U, 
no further precision in specifying the intermediate state's energy is needed, and the leading 
order term is sufficient. When (U — u) becomes of order t (resonant regime), all of the terms 
in the expansion in t/(U — u) are of the same order, and the leading term can no longer 
provide a solution. However, it is in this regime that the cross-secions sensitively depend 
not just on the magnetic, but also on the carrier properties, and this makes it a especially 
important to the field of high-T c superconductivity. 

We will discuss the form of the Raman vertex in both the nonresonant and resonant 
regimes later in the paper and here merely note that, on physical grounds, the multiple hops 
affect the single-particle Green's function in two distinct ways. First, the multiple hops is 
a way to generate the band structures of the doubly occupied sites and holes, i.e., they 
produce a dispersion in the coherent part of the single-particle Green's function. Second, 
they produce a distortion of the spin background around each hole which gives rise to the 



incoherent part of the Green's function. Since resonance is a phenomenon contingent upon 
the vanishing of the denominators in perturbation theory, it is plausible to assume that 
the incoherent parts of the Green's functions are of lesser importance. We thus model our 
particle and hole Green's functions as fully coherent, but having a certain dispersion, which 
at half-filling we take to be that of the spin density wave (SDW) solution of the Hubbard 
model.ii We regard our ability to address the experimental results on the basis of this 
dispersion relation as partial evidence of the validity of the SDW picture not only for the 
antiferromagnetism, but also for the hole dynamics. 

We wish to stress that this is not a trivial distinction. There have been many interest- 
ing theoretical studies of the single hole problem in the Hubbard and the t-J models.@~@ 
However, even such a basic feature as the band structure of a single hole has not been 
experimentally tested. Indeed, in the low-doping regime the holes are trapped near the 
dopant sites. Even in the translationally invariant case, the phonon polaron effect JH 
could still destroy the hole dispersion predicted on the basis of purely electronic models. 
Resonant Raman scattering right at half-filling is sensitive to the details of the dispersion of 
the fermionic quasiparticles and thus provides an experimental information about the single 
hole dispersion without having to dope the material. 

Before proceeding to the description of our calculations, we would like to make two 
comments about the model we use. First, in this paper we exclusively consider the one-band 
Hubbard model. This by itself is a simplification, since in the "first-principles" calculation 
one would start with a three-band model for the C11O2 unit. Within this model, the lowest 
gap (about 2 eV) is a charge-transfer gap between Cu and O bands.S However, it is generally 
accepted that at energies comparable to this gap the hybridization between Cu and O 
orbitals is relevant, and one can effectively describe the system by a single degree of freedom 
per CuOi unit, which in turn implies that the three-band Hubbard model can be reduced 
to an effective one-band Hubbard modeli^S The relevant parameter which allows the 
reduction in the number of degrees of freedom is the splitting of the triplet and singlet states 
of the O and Cu spins in the Zhang-Rice theory.0 Their estimate of this splitting is several 



electron volts, and thus the singlet states relevant to our analysis are well-separated from 
the triplet states. This is corroborated by the measured behavior of the dielectric constant. 
For example, in GdCuC>2 there are strong absorption features in the imaginary part of the 
dielectric constant around 1.7 — 2.7 eV and around 5 eV, but only a featureless continuum 
in between. We therefore expect that the one-band Hubbard model already captures all 
the essential physics of Raman scattering up to photon energies of order 3 — 4 eV. 

Our second comment concerns the hopping term in the Hamiltonian. In this paper, 
we restrict our model to contain only the nearest-neighbor hopping, t. Meanwhile, it was 
suggested!! on the basis of the photoemission data for YBaCuO-j, that the next-nearest- 
neighbor hopping term in the 123 compounds is rather large, t' ~ —0.51 For La-based 
compounds, t' was estimated!! to be t' ~ — 0.2t. At finite doping, the restriction to only 
the t term requires care, as the spin-density wave theory (which we will use in most of 
our considerations) is meaningless near the minima of the quasiparticle band, unless one 
takes into account the self-energy corrections which lift the degeneracy of the quasiparticle 
spectrum along the reduced Brillouin zone boundary.il'0 Fortunately, as we will see below, 
the two-magnon Raman scattering in cuprates is not a phenomenon that is associated with 
the minima of the quasiparticle band, and therefore the details of the quasiparticle dispersion 
right near those minima will not be important in our analysis. For this reason, we believe 
that for a qualitative analysis of Raman scattering in cuprates, one can use the simplest 
model with just the nearest-neighbor hopping. 

The paper is organized as follows. Sec. II discusses the general formalism of the two- 
magnon Raman scattering in antiferromagnetic insulators. In Sec. Ill, we consider the 
nonresonant scattering. In particular, we show how the results of Shastry and Shraiman 
are reproduced in the SDW formalism. We will also discuss the location and shape of the 
two-magnon peak. Sec. IV is the main part of the paper. We argue that the phenomenon 
of multiple resonance, described in detail in that section, is a way to explain the unusual 
features in Figs. [TJ and 0. We study the singularities of the integrals corresponding to 
that type of resonant behavior, and present the relevant calculations. Finally, in Sec. V, 
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we present a comparison with the experimental results and state our conclusions. Some 
technical details of our calculations are contained in the Appendix. Some of the results of 
this work have been already presented in a short paper .0 



II. GENERAL CONSIDERATIONS 

We start with the general theory of magnetic Raman scattering in the Hubbard model. 
The one-band Hubbard Hamiltonian is 

H= -t ( c L c i,- + h.c.) + U^ni^. (2) 

<i,j> i 

A straightforward procedure to derive the coupling of light to the fermions in this model was 
recently described by Shastry and Shraiman.El In the presence of the slowly varying vector 
potential A(x, t), each fermion operator acquires a phase e ( ie A c )J A ' dI ) anc i the hopping term 
gets transformed into 

fft = -tE A ^V,,0, • h.c. (3) 

<i,j> 

The vector potential is supposed to vary slowly on lattice scales, and one can then approxi- 
mate the phase factor by 

J' A(l) • rfl « A ■ R^, (4) 

Substituting this into the expression for H t and expanding in powers of the vector potential, 
one obtains upon transforming to the momentum space 

Ht = Ht=° - £ E h • A- fl + \{^-f E E A- 91 • W ■ A_ 92 , (5) 

fit- q flC qi q2 

where the current and the stress-tensor operators are 

■ a _ sr^ t 

3q ~ QJ, C k+q/2,a C k~q/2,a, 

k a 
52, 



a,/3 _ \ "* ® e fc t (a\ 



and ek is the electron dispersion, 

ek = — 2t(cos k x + cos ky). (7) 
The vector potential can be quantized in the usual way 

A 9>A = (/ 9 (e A Z_, + 3$), (8) 

where g q = (2%hc 2 / (cUgV)) 1 ^ 2 and cu q is the light frequency. The velocity of light is several 
orders of magnitude larger than the Fermi velocity, and therefore we can safely set the 
momenta of photons equal to zero. 

Consider now separately the last two terms in Eq. (|5|), which contain the vector potential. 
The last term is already quadratic in A, so for photon scattering we have to consider only the 
first-order contribution from it. This contribution does not lead to resonant scattering and, 
at least in the resonant regime, can be neglected compared to the resonant scattering from 
the first term. For this reason we neglect the last term in Eq. ([5]) in this paper, and focus 
only on the first term which leads to photon scattering in the second order of perturbation 
theory. This second-order scattering process involves an intermediate particle-hole state of 
the fermionic system which can emit or absorb collective bosonic excitations before collapsing 
into the outgoing photon. In particular, in the two-magnon resonant scattering, a photon 
with the energy u>i is injected into the antiferromagnetically ordered set of electrons. This 
photon creates a virtual particle-hole pair which then emits two spin-waves with momenta 
k and — k, and then annihilates into an outgoing photon with the energy ujf. 

Ignoring the polaritonic effects, the Raman scattering cross sectionEl is obtained from 
the lowest-order Golden Rule 

87r 3 e 4 

R = TMF2 E I < f\ M *\ l > I' 5 ^ ~ huJ f + e * - e /)> ( 9 ) 

ri V ujjUj j ^ j 

where e« and ej are the energies of the initial and final states of the system (ej — 6j = Q 
is a total frequency of two magnons in the final state), and the summation over the final 
states / includes the integration over the photon momenta. Further, < /|Mr|z > is a matrix 
element, given by 
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M f\M x-\ < f\)kf-e*j\ n > <n\i- ki -e i \i> < ■e l \n> < n\j kf ■ e} \i > - 

M R =< f\M R \l >= — Z -TT 1 Z —? 

n l €i + hUi - e n + id €i - e n - hujf + to 

(10) 

Here the summation is over the intermediate electronic states, labeled as n. Our primary 
goal will be to calculate the dependence of this matrix element on the incident photon 
frequency. 

For completeness, we also list a number of possible experimental scattering geometries.il 
They differ in the polarizations of the incident (e$) and scattered (e/) photons. For the 
linearly polarized light, the scattering geometries are 



Big 





x + y 




x + y 




y/2 ' 


g / = 


V2 




x + y 




x-y 


e, 


V2 ' 


e / = 




e, 


= x, e f 


= y- 





B 2g : e { = x, e/ = y. (11) 

For the circularly polarized light, the scattering geometries are 

x + zy „ x + zy 

rD . x + zy A x-zy 

M:e i = - ;rj e, = -^-. (12) 

The scattering in each geometry measures a particular combination of the components of 
the scattering tensor. It is convenient to decompose < f\M R \i > into four one-dimensional 
representations of the 2D square lattice symmetry group : 

= (x\M R \x) + (y\M R \y), 
= (x\M R \y) - (y\M R \x), 
Ml 1 = (x\M R \x) - (y\M R \y), 

= (x\M R \y) + (y\M R \x). (13) 

Then, the scattering in the experimental A\ g geometry measures a combination of M^ 1 and 
M^ 2 , the scattering in the Bi g geometry measures a combination of M^ 1 and M^ 2 , etc. In 
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the two-magnon Raman scattering, the energy transfer Ui — uoj is small compared to the gap 
between conduction and valence bands. In this situation, the dominant components of the 
scattering tensor are found experimentally to be those of M^ 1 and M^ 1 symmetry^! It is 
therefore sufficient to study scattering in the Ai g and B\ g experimental geometries. 

The states in the general expression of Eq. (jT^) for the matrix element are the full 
many-body states of the system, which contain both the electronic and spin excitations. 
Therefore, to fully describe the two-magnon Raman scattering we will need two different 
types of vertices. The first is the interaction between the vector potential of light and the 
fermionic current density - this vertex is given in Eq. fl5|). The second is the interaction 
between the fermions and the spin-waves. 

A convenient way to express both types of vertices on equal footing is to use the 
spin density wave (SDW) formalism!! to describe the electronic state at half-filling and 
the excitations around it. In the SDW formalism, one introduces a long-range order in 
S g = J2k ^k+q^oPa^Ck^ with q = Q = (tt , 7r) and uses it to decouple the Hubbard interac- 
tion term. The diagonalization then yields two bands of electronic states (the conduction 
and valence bands) with the gap 2A ~ U in the strong coupling limit that will be assumed 
throughout this work. 

In terms of the conduction and valence band quasiparticle operators a\ a and b\ a , the 
quadratic part of the Hubbard Hamiltonian takes the form 

H = Y^E h {a^a ha -b\ a b 1m ), (14) 

ka 

where the prime restricts summation to the magnetic Brillouine zone, and the quasiparticle 
energy is E k = \J e\ + A 2 . We will also need the relation between the new and the old 
quasiparticle operators. It reads 

Ck,cr = U k Cl k)Cr + Vkbk }C r, 

Ck+Q,a = sgn(a)(u k b kj<7 - v k a k:<7 ), (15) 

where the fermionic Bogolyubov coefficients are 
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Note that in the large-U limit, 2ukVk — 1. 

We now discuss the vertices. The current operator can be rewritten as 



h=0 — 2^ ~Qfc C kcr C kcr ~ 2^ \ C kcr C kcr ~ c k+Q,a C k+Q ,cr ) ■ I 1 ' J 



Upon performing the Bogolyubov transformation, we obtain the current in terms of valence 
and conduction band fermions, 

j g =o = Y!^[( 2u k v k)( a l* b ka + bl a a ka ) 

kcr 

+ ( U l - V k)( a L a ka - &!A<t)]- ( 18 ) 

Note that the terms with only valence or conduction fermions have a smallness mt/U from 
the Bogolyubov coefficients. Keeping only the leading terms, we obtain 

j g =o => Y!-gr( a L b ka + b[ a a ka ), (19) 

kcr 

which is the expression we use almost exclusively in this paper. 

We can now obtain an expression for the optical conductivity <r(q = 0,u) in the SDW 
state. From Kubo formula, in terms of the particle current given above, 

a(u) = e 2 -J2S(e n -e Q -u) \(n\j x (q = 0)|0)| 2 , (20) 



where e n — eo is the excitation energy between the ground and excited SDW states. At 
half-filling, we only have to retain the intraband terms in Eq. flT8p, and, using also Eq. (|TJ|) 
for the valence and conduction fermion energies, we obtain 

^)=e 2 — ^'sin 2 ^ ^ S(u-2E k ), (21) 

U ka ^k 



which agrees with the result in Ref. |23j. We have obtained a compact closed form answer for 
<j(uj) in terms of complete elliptic integrals. The answer is nonzero when the photon energy 
is within the bands, 
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A 



4 



yf\K J 1 



A 
4 



) 



< A < 4, 



(22) 



where A = (u 2 — (2A) 2 )/16t 2 counts the photon energy from the optical gap edge, and 
we momentarily restore fi for reference purposes. When t/U <C 1, 2 A/a; m 1, and A ~ 
(u — U)U/8t 2 . The plot of a{ui) at t -C U is given below in Fig. |ll](a). It has a square root 
singularity at A = + and vanishes as A — > 4~. The singularity at A = + is due to the 
degeneracy of the mean-field dispersion Ey. along the reduced Brilliouine zone boundary. As 
we described above, this singularity is an artifact of mean-field SDW approach. When the 
degeneracy in Ey. is broken by either quantum fluctuations or hopping to further neighbors, 
only the more usual logarithmic singularity will survive at uo = 2 A. On the other hand, 
the vanishing of ct(uj) at A = 4 is due to the vanishing of the current form-factor sink x in 
Eq. (|2~ID at k = 0, which is where the maximum optical gap occurs for SDW states. There is 
no mean-field degeneracy at k = 0, and thus o = at the top of the fermionic band should 
survive also beyond the mean-field approximation. 

The remaining ingredients for Raman scattering are the magnon propagator and the 
magnon-fermion scattering vertex. In the SDW formalism the magnons are described as 
collective modes in the transverse spin channel.0 Specifically, the spin-wave excitations 
correspond to the poles of the transverse spin susceptibility 



The total transverse susceptibility is given in the SDW theory by the diagrams containing 
a sequence of bubbles made of the conduction and valence fermions with the four-fermion 
vertices connecting them. The restriction to only such diagrams can be justified if one 
extends the original S = 1/2 Hubbard model to large S by considering 2S orbitals at a 
given site@ At half-filling this extension transforms the Hubbard interaction term into the 
Hund's rule coupling which favors the maximum possible spin S at each site. Notice that 
in this situation the mean-field gap between lower and upper Hubbard bands is related to 
U as 2A = 2US. In relating the large-S" calculations to the S — 1/2 case, it is the gap A 




(23) 
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which should be kept fixed. 

The SDW spin susceptibility has been considered several times in the literature. HElIll 
Since the unit cell is doubled due to the presense of the antiferromagnetic long range order, 
we have two susceptibilities — one with zero transferred momentum, and one with the 
momentum transfer Q = (ir,7r). The explicit forms of these susceptibilities are 

1 



-s 



X ± (q,q + Q,uj) = -S 



' l-7g 1 L 

1 + 7 9 UJ — VL q + i5 UJ + VL q — i5 

1 1 



to — Clg + iS u + Cl q 



iS 



(24) 



Here fl q = 4JSyl — 7^ is the magnon frequency, and J = 4t 2 / (2S) 2 U is the exchange 
integral. 

Further, a sequence of bubble diagrams can be viewed as an effective interaction between 
two fermions mediated by the exchange of a spin- wave. The spin- wave propagators are 
i(Te g (t)e^(0)) £J — (Q q -oj - id)' 1 and i(Tejj(t)e g (0)) w = (fl, + u - z<5)~\ where e q (e q ) are 
the boson creation (annihilation) operators, subindex to implies Fourier transform, and the 
momentum q runs over the whole Brillouin zone. A simple experimentation then shows that 
the forms of the two susceptibilities are reproduced if one chooses the following Hamiltonian 
for the interaction between the original fermionic operators and the magnons: 

2A 

k q 



H, 



el— mag 



[Vq C k+q,a C k,l3( e -q + e s) + Vqcl+q^k+QA^-q ~ e q) 



5, 



a ,-p- 



(25) 



Here and below a prime indicates that the summation is over the reduced Brillouin zone 
The expressions for r] q and TJ q are 

Performing now the Bogolyubov transformation, we obtain the Hamiltonian for the in- 
teraction between the magnons and the conduction and valence band fermions 
/ 

H e l- m ag = ( a lak a il3,k+ q e\ $aa(^, ?) + Aak h ip,k+q & \ ®bb(k, q) 

k q 

+ a Lk b if3,k+ q el 3> ab (k, q) + bl ak aip ik+q el <$> ba {k,q) + H.c. ) 6 a> ^. 
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(27) 



To leading order in t/U, the vertex functions are given by 

1 



$aa,w>0, q) = ±(efc + e k+q )r] q + (e k - e k+q )r] 



2S 



®ab,ba(k, q) = 2A [ Vq =F r/J -L. (28) 



We are now in a position to proceed systematically with the diagrammatic formulation 
of the two-magnon Raman scattering. 

III. NONRESONANT SCATTERING 

A. Raman matrix element 

We start with the situation when the photon frequencies are much smaller than the gap, 
A ~ SU. In this case, it has been been shownll that the Loudon-Fleury Hamiltionian0 gives 
a proper description of magnetic Raman scattering. We first show how this Hamiltonian 
can be reproduced in our momentum-space diagrammatic formalism. The key points in the 
derivation are the following: (i) to leading order in t/U the current operator necessarily 
transforms a valence fermion into the conduction one and vice versa; (ii) the vertex strength 
for a magnon emission accompanied by the fermion scattering from the valence to conduction 
band (and vice versa) is of order U, while for the fermion scattering within either the valence 



or conduction band, the vertex is of order t, i.e., much smaller (see Eq. (28) ); (iii) for 
incoming/outgoing photon frequencies smaller than the gap, all the denominators in the 
diagrams are of the order of U. We can then identify three simple diagrams contributing 
towards the matrix element Mr of Eq. fllPf) to leading order in t/U. They are shown in 
Fig. [|. Notice that while the diagrams of Figs. [|(a),(b) have a superficial resemblance to 
those encountered in the theory of two-phonon Raman scattering in semiconductors]^ the 
diagram of Fig. |](c) looks rather different. We will see shortly that the extra diagram is 
needed for the answer to have correct symmetry properties. 

Performing the internal frequency integration in these three diagrams, we obtain 
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•vL l, = -4/, A E'^-^W^i-e 



4A 8A(4A 2 + Q 2 



L R 



dk 



dk f J \AA 2 -tt 2 (4A 2 - tt 2 Y 



(29) 



for the diagram in Fig. |](b), 



'dt k 
dk 



'de k \8A(4A 2 + ft 2 



9k f ) (4A 2 -n 2 y 



(30) 



for the diagram in Fig. |](c), and 



fofc- g \ / 4A 
9k ' e/ j V4A 2 -fi 2 



(31) 



for the diagram in Fig. |](a). In these three expressions Q is a frequency equal to Ui or u;^. 
Specifying it more precisely at this stage is devoid of meaning since we would then have to 
take into account subleading diagrams in the t/U expansion. We also defined 



Vq 



2 \ . I\ - 7 2 



+ 1 



\lq\ 



1"7 2 



These are related to r) q and r\ via 



^Vq = Vq + Vq] ^\ = (Vq ~ Vq) ■ 



(32) 



(33) 



Adding up the above expressions, and taking into account that 



?bk- e ^ 



d_e k=q 
. dk f , 



= J2'(- 2t ) 2 [(xsin/c x + y smk y ) ■ 6*] (±sm(k x - q x ) + y&m(k y - q y )) ■ ej 
= t 2 [e^e^ cos q x + e^e^ cos ^ 

we obtain for the Raman matrix element 
2A 



M. 



R 



4A 2 - u 2 



(34) 



Six 6 fa; + e iy e* f )2fi q X q - (e 



(35) 



We now demonstrate that this equation coincides with the Loudon-Fleury vertex. For 
this, we first rewrite the Loudon-Fleury Hamiltonian of Eq. (|I|) in terms of Holstein- 
Primakoff operators to leading order in 1/S 
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r,fi={x,y} 

a 



77 [^/^CirrS,- ' (^r+x + S r _f) + e^ejyS r • (S r +,, 



fy^y 

cos q x (a q a_ q + a)a}_ q ) + 2a£ 



= 2Sa Y!{ e }x e i 
g 

+ e* fy e iy [cos q y (a q a_ q + a\a ] _ q ) + 2a\a q } 



q a Q 



(36) 



We then apply the Bogolyubov transformation to the actual magnon operators that diag- 
onalize the Heisenberg model to first order in 1/S. This transformation involves the same 
factors X q and fi q that we had introduced earlier (in Eq. fl32|)) and has the form 



f^qO — q 



(37) 



At T = the only relevant terms in the Loudon-Fleury Hamiltonian are those containing 
two magnon creation operators, and restricting ourselves to those, we obtain the matrix 
element from the Loudon-Fleury Hamiltonian in the form 



2<\S{< f r ( , 



+ e fy e iy 



cos q x (\ 2 q + fi 2 q ) -2fx q \ q 
cos q y {X 2 q + fi 2 q ) - 2/^A ? ] } 



(38) 



Comparing the two expressions for the Raman vertex, Eqs. (|3^) and (^), we observe that 
they coincide if one chooses a = 16t 2 A/[2S(AA 2 — u> 2 )]. Apart from the factor 25* (= 1 
for S = 1/2), this is exactly the expression which Shastry and Shraiman obtained in their 
derivation of the Loudon-Fleury vertex for the Hubbard model.0 

Further, an examination of our diagrammatic derivation shows that the terms in Eq. fl35|) 
containing cos q x and cos q y correspond to the S + S~ terms in the Loudon-Fleury Hamltonian, 
while the rest of the terms come from the S Z S Z part. Going back to Eqs. (|29|), fl3"0|), and 
(pl|) and the diagrams they correspond to, we can see that the S + S~ terms come from the 
diagram of Fig. f|(a), which emits one magnon from the hole line and other from the electron 
line. The S Z S Z terms come from the diagrams of Figs. |](b),(c), which emit both magnons 
either from the electron or the hole side of the fermionic bubble. Notice that the "unusual" 
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diagram of Fig. |](c) canceled one of the terms in the one of Fig. f|(b). Without that 
cancellation we would not have obtained the correct Loudon-Fleury expression for the S Z S Z 
term. This in particular would have resulted in a violation of the selection rule prohibiting 
the Ai g scattering. 

There is an interesting parallel between this cancellation and one that occurs in the 
theory of two-phonon Raman scattering in semiconductors.0 In that case, there is a large 
cancellation between the two-phonon terms obtained by the iteration of a single-phonon — 
electron vertex, and the diagrams containing a single two-phonon — electron vertex. The 
cancellation in the phonon case is a consequence of the translational symmetry. In our case, 
the diagram of Fig. (|(c) can be seen as effectively containing a single two-magnon — fermion 
vertex. The partial cancellation due to this diagram restores the spin rotation invariance 
of the Loudon-Fleury Hamiltonian. It is then reasonable to assume that some cancellations 
may occur in higher orders in t/U, though we did not check this in explicit calculations. 

Despite the agreement obtained between the two approaches, we are not done yet. We 
have so far been performing an expansion in t/U, and there exist other graphs which are 
of the same order in t/U as the ones retained. Examples are given in Fig. [5|. However, it 
is not difficult to observe that all these extra diagrams have a smallness in 1/S. Indeed, 
each fermion- magnon vertex has an overall factor of 1/V2S (see Eq. (p8|)). Associated with 
each internal magnon line, we have two such vertices, and no additional summation over 
the orbitals in the internal fermionic loop. Altogether, inclusion of an internal magnon line 
results in an extra factor of 1/25*. This completes our demonstration of the equivalence of 
our formalism to that of Shastry and Shraiman in the non-resonant region. 

B. Two-magnon peak 

In the preceding subsection, we calculated the Raman matrix element Mr for the emission 
of two magnons. Here we consider the location and shape of the two-magnon peak. We first 
observe from Eq. (^) that within the Loudon-Fleury picture the scattering in the A lg 
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geometry (e^ = ^jf, &f — vanishes because the Raman vertex M R is equal to zero. 
This result can be expected because at half-filling (and at large U) the Hubbard model 
is equivalent to the Heisenberg model, and in the A\ g geometry the Loudon-Fleury and 
Heisenberg Hamiltonians commute with each other. On the contrary, in the B\ g scattering 
geometry (ej = e/ = ^jtf) the Raman vertex, as obtained from Eq. (^), is finite 

Mf 1 ~ + A^cosg, - cosg,) = (39 ) 

v 1 - 7 * 

Substituting this into the expression for Raman intensity, Eq. (||), we obtain 

RBt (w) „ ^ (eos, -eos,,)^ _ 2SJ?) _ (4Q) 

where u = Ui — ujf. This expression has a divergence near uj = 8JS, since the magnon 
spectrum is flat and the density of states infinite at the Brillouine zone boundary. It is 
well-known, however, that this result changes qualitatively when one takes into account 
the effect of the magnon-magnon interactions in the final state. This was first explained 
by Elliott et. a/.0 and experimentally verified by FleuryS for RbMnF3, which has a cubic 
structure. The calculations were further extended by Elliott and Thorpe.@ They were per- 
formed for the two-dimensional case by Parkinson,^ and then verified for K 2 NiF 4 by Fleury 
and Guggenheim.0 Since a rather up-to-date treatment of the role of the magnon-magnon 
interactions is available in the literature,i@ we merely sketch the relevant calculations with 
an emphasis on the 1/5* expansion, which treats the magnon-magnon interactions in a sys- 
tematic way and which, to our knowledge, has not been discussed before in this context. 

The standard way to treat the effects of multiple magnon-magnon scattering in the 
Raman problem is to write the Heisenberg Hamiltonian in terms of the Holstein-Primakoff 
bosons, diagonalize it, and keep only the interaction term ~ ar k /3_ k /3-iati, which is responsible 
for multiple scattering of two magnons in the vacuum (Fig. ^]). This is the procedure which 
Parkinson^ and others used to derive the two-magnon Raman intensity. We note however 
that for arbitrary S, the restriction to only a single interaction term is not justified, as after 
the diagonalization of the quartic term in the Holstein-Primakoff bosons, one also obtains 

19 



processes of the form, say, a^-fc/^ljCKj , which also contribute to the Raman intensity. The 
way to avoid this complication is, again, to study the large- S limit. The key point is that 
at infinite S the Raman intensity is divergent at uj = 2Q q = 8JS, which is the maximum 
possible two-magnon frequency. It is achieved when both magnons are right at the Brillouin 
zone boundary. At this boundary 7^ = (cosk x + cosfc J/ )/2 = 0, and the anomalous term in 
the bosonic quadratic form in the Heisenberg antiferromagnet vanishes. In other words, the 
antiferromagnetic magnons right at the Brillouin zone boundary behave as free particles. 
Now, if S is large, the shift in the two-magnon Raman peak position due to the magnon- 
magnon interactions is small, so that one can still consider magnons as free particles and 
neglect the anomalous part in the quadratic form. In this situation, the only interaction 
term has two creation and two annihilation operators. Moreover, in the Raman problem the 
total momentum of two magnons is zero, and we can therefore write the interaction term in 
the truncated form 

H^t = -^E'E'Tk-jat/^-jaj- ( 41 ) 

k I 

The prime stands for summation over the magnetic Brillouin zone. We can further decom- 
pose 

lk-i = Ikli + Ikli + + %% (42) 
where the different symmetry factors are 

Ik = ^ (cos k x + cos ky) , 7 fc = i (cos k x - cos k y ), 

1 1 
l k = -(sin k x + sin k y ), 7 fc = -(sin k x - sin k y ). (43) 

For the B ig scattering geometry, Mr oc j q , and hence only 7^77 survives at each vertex. The 
remaining calculation of the series is straightforward, and we obtain for the Raman intensity 

I 



where 



1 + 4^ 



(44) 
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_ J_ w (cosg x -cosg^) 2 

and we introduced the reduced frequencies uj = u/SJS and fl q = flg/AJS. 

The imaginary part of Eq. (f45[) can be expressed in terms of complete elliptic integrals, 



Im I 



7T\/1 — UJ 2 



[2-uj 2 )K(uj) -2E(Q) 



(46) 



One can then express / via a dispersion relation, 



8 ^,,[(2-,^M-2 gM ] 
V ; Wo v 7 ! 3 ^ (a) - v + i<?) V ; 

This one-dimensional integral is particularly convenient for a numerical evaluation of /. 

For S = oo the denominator in Eq. ([0]) is irrelevant, and this expression reduces to 
(fjQ|), with the only difference in that we have set 7 g — > in the Mr of Eq. fl39|), in order 
to be consistent with the approximations made in the derivation of the magnon-magnon 
interactions. 

For large but finite S, the peak in the Raman intensity is located close to but not exactly 
at uj = 1. Expanding the integral in Eq. fl45|) in 1 — u, we obtain I = A + iB where 



l-VT 



-,2 



A = -4 , B = 7 _= log(l - O 2 ). (48) 

V 1 — oj n v 1 — v 

Substituting this into Eq. (|44] ) we obtain after simple manipulations 



VT^a? iog(i-cu 2 ) 

\0g 2 (l -CU 2 )+7T 2 (l-(5+l)(l - U)f ' 



iW^oc ^ i;r ,L (49) 



We see that when the magnon-magnon interaction is included, Rb 19 (^) n o longer diverges but 
rather has a maximum at \/l — u 2 rs 2 log (S 1 + l)/(n(S + 1)). Moreover, right at the zone 
boundary Raman intensity turns to zero. If we formally set S = 1/2 in our 1/5* expression 



for Raman intensity, Eq. fl49|) , we obtain a peak centered at u — 0.696 (or to = 2.78 J). 



The peak position is close to the uj = 2.92J result obtained by Canali and Girvin.B They 
used the same expression for I as we did and restricted to the same sequence of bubble 
diagrams, but they did not make a formal assumption of large S in their calculations of 
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intensity and therefore did not expand I in 1 — u and set 7 g = 1 in fl40|). The agreement 
between our results and theirs shows that the two-magnon peak position is rather robust. 
On the other hand, the profile of the two-magnon peak depends strongly on at which stage 
in the calculations one extends a formal 1/5* expansion to S — 1/2. We found that the use 
of the 1/S expansion results for A and B yield much broader peak than one obtains using 
the exact expression for I. This last form of the two-magnon profile is consistent with the 
results of numerical calculations^ and we therefore believe that it is closer to reality. The 
profile of the two-magnon intensity is plotted in Fig. 0. 

Canali and Girvin also considered the effect of including the spin-wave velocity renor- 
malization factor into the magnon propagator. To first order in 1/S, this renormalization 
factor effectively shifts the value of J to J e ff = J(l + 0.16/2S 1 ). The same renormalization 
must indeed be present in our approach, so the position of the peak in our theory is at 



The peculiar features of the two-magnon profile were first reported by Elliot and Thorpe-^ 
and Parkinson,^ who performed the RPA summation without referring to the 1/S" expan- 
sion, neglected spin-wave velocity renormalization and also kept the interaction terms with 
four creation and annihilation operators, which were assumed to have the same overall factor 
as in Eq. (|4l|). For S = 1/2, a numerical solution of their equations yields a narrow peak 
at to = 0.675 (i.e., uj = 2.7 J), which is not far from both our result and that of Canali and 
Girvin, though we believe that the restriction to a single interaction term is better justified 
as long as the virtual magnons are located near the Brillouin zone boundary. 

Experimentally, the two-magnon peak has been observed in Lg^CwC^ at 3000 cm -1 (see 
Ref. §), which for uj ~ 2.78J e // ~ 3.22J yields J ~ 0.116 eV which is roughly consistent with 
J ~ 0.125 eV as inferred from the neutron scattering^ and NMR data.@ The Parkinson 
result corrected by the spin-wave velocity renormalization (uj = 2.7J e ff ~ 3.13J) gives 
somewhat better agreement with the neutron data, but the difference between the two 
values of u is rather small and practically irrelevant to subsequent analysis. In the rest of 
the paper, we will simply refer to 2.8J e // for the two-magnon peak position. 



uj ~ 3.22J. 
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We now proceed to the central topic of our paper, which is the discussion of Raman 
scattering in the resonant regime. 



IV. RESONANT SCATTERING 



A. Raman matrix element 



In the section on non-resonant scattering we proceeded systematically in filtering out 
the diagrams that were small in t/U or 1/S. Our very first step was to observe that the 
magnon-fermion vertex was of order U when the fermion scattered between the valence and 
conduction bands, and of order t when it scattered within either the valence or conduc- 
tion band (see Eq. (p8|)). As the terms in the denominators obtained upon the internal 
frequency integration were always of order U, all the diagrams with intraband scattering at 
any magnon-fermion vertex were small by powers of t/U. 

It turns out that this property of the denominators is no longer true in the resonant 
region, where the photon frequencies u>i and Uf differ from the gap 2 A by quantities of 
order JS. As a result, some previously omitted diagrams become, as we will see, not only 
important, but actually dominant. 

Let us first consider the form of the Raman spectrum without the final state interactions. 
Then we have 



Consider now the graphs with only the intraband scattering at the fermion-magnon vertices, 
Fig. || With all the multiplicity factors included, the internal frequency integration results 
in 



R(u) oc \M R \ 2 S(ui -Uf- 2Q q ) 



(50) 






(51) 



R — 



(ui - 2E k + i5)(ui - Q q - E k - E k _ q + i5)(u f - 2E k _ q + iS) 



for the graph of Fig. ||(a), and 
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M fl - op. ,,-j;v n _ P. _ P. i.-AV OP. x,^ ^ 



k 



[UJ 



2E k + i5){ui -n q -E k - E k „ q + iS)(co f - 2E k + iS) 



for the graph of Fig. §(b). The relation Ui — 2Q q = ujf is to be remembered here. 

Without performing the integrals, the naive order of magnitude estimate is obtained in 
the resonant region by assuming that the resonant terms in the denominators are of the order 
of the bandwidth, i.e., O(JS). Then the two diagrams above are of order t A /(JS) 3 . The 
diagrams that were dominant in the nonresonant region have fewer resonant denominators. 
Going back to Eq. fl35p for the previously considered graphs, we find that they are of order 
t 2 /(2A — uj) ~ t 2 /JS. We see that in the resonant regime, the new diagrams are larger 
by a factor (t/JS) 2 ~ U/J. However, this estimate, while correct for the large-S* case, can 
be somewhat misleading at the physically relevant value of S = 1/2, as it does not take 
into account strong self-energy and vertex corrections which are relevant at u ~ 2A, even 
though we are considering the situation exactly at half-filling. For example, the leading order 
vertex correction to the magnon-fermion vertex, shown in Fig. |^(a), is of order t(U/JS) 2 , 
whereas the bare vertex is of order t. Simultaneously, the self-energy correction in Fig. [|(b) 
contributes an extra factor of (U/JS). Both corrections are small only if we require that 
U/JS <^ 1, which we indeed do not expect to be satisfied for S — 1/2. 

The form of the low-energy theory at U/JS ^> 1, has been discussed by a number of 
authors. For our considerations, it is essential that both the numericali! and variationafi 
studies have found that the resulting quasiparticle Green's function G(k, u) has a quasipar- 
ticle pole at low energies, and that the fc-space dispersion of that quasiparticle has a width 
of order JS. Self-consistent calculationsSfill! also demonstrated that the product of the 
quasiparticle residue Z and the effective (renormalized) interaction between quasiparticles 
U e ff scales as J (without corrections this product was equal to Ux 1 = U). In this situation, 
the relative factor between the new and the old diagrams is U e ffZ/J = 0(1). The order 
of magnitude estimates therefore only point to the new diagrams as being comparable in 
importance at resonance to those that fully dominated away from the resonance. 

Fortunately, this naive order of magnitude estimate is not the whole story. When the 
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integrals are actually performed, they may (and as we will see, they do) yield singular answers 
for some photon and/or magnon frequencies. In that case, the most singular integral will 
dominate. Notice that because the singularity comes only from the coherent part of the 
quasiparticle spectral weight A(k, uj), we do not need to consider incoherent part of A(k, lu), 
which by itself can be substantial and spread over the energy scale exceeding J (some studies 
predict that A(k,u) spreads over the scale of the order J (j) 7 , where 7 is a number < 1). 

We now consider in detail the singular behavior of various diagrams. The diagrams which 
contribute to the Loudon-Fleury Hamiltonian have only singly resonant denominators that 
lead to singularity in Mr when either the incoming or the outgoing photon frequency coin- 
cides either with the gap, 2A, or with the top of the band.0 Near the gap, the integration over 
the intermediate momenta yields Mr ~ \uJij — 2A| -1 / 2 for the mean- field form of the quasi- 
particle dispersion, or an even weaker logarithmic singularity Mr ~ log — 2A| for the 
renormalized quasiparticle spectrum, which has minima at (±7r/2, ±7r/2) and a parabolic 
dispersion around the minima. This last form of the fermionic spectrum was obtained in the 
numerical, variational, and perturbative studies of the Hubbard model beyond the mean- 
field level.ilSB Near the top of the band (k rs 0) , the would be logarithmic singularity due 
to the vanishing of the denominator is in fact absent because ^ in the numerator of Mr 
also vanishes at k = 0. 

Consider next the diagram in Fig. |8](a). Its analytical expression is given by Eq. (|51~1). 
The denominator in Eq. (|5T|) has three terms, one of which is half of the sum of the other 

twoizl 

uj l -Vt q -E k - E k+q = - [(ui - 2E k ) + (ujf - 2E k+q )} . (53) 

Therefore, if we find k and q such that u>i = 2E k and ujf = 2E k+q , all three of the denom- 
inators in ([[1]) vanish simultaneously. This is known as triple resonance, which is a partic- 
ular case of the more general notion of multiple resonance known in Raman scatteringliE 
Clearly, a necessary condition for a triple resonance is that the fermionic bandwidth be larger 
than the magnon bandwidth. This condition is satisfied in the mean-field theory, where the 
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magnon bandwidth is 4JS, while the fermionic bandwidth, found from 

E k ~ A + 2JS(cos k x + cos k y ) 2 , (54) 

is equal to 8JS. The bandwidths of "dressed" Fermi and Bose quasiparticles are indeed 
somewhat different, but we assume that the fermionic bandwidth is still larger than the 
magnon bandwidth. 

We will present a detailed study of the momentum integration near the triple resonance 
in the next subsection, and here merely note that, unlike its single counterpart, the triple 
resonance is not tied to the minimum of either E^ or Ej, +q . Thus, unlike for the Loudon- 
Fleury terms, for triple resonance the difference between the mean-field and the renormalized 
forms of the quasiparticle spectrum does not yield qualitatively different answers for the 
singularities. For this reason, we will keep working with the mean-field form of E k . In the 
next subsection we will find that triple resonance does indeed yield the strongest divergence 
of the Raman scattering cross-section. 

Finally, in the diagram of Fig. |](b), we have a product of (u>i — 2E] C ), (u>f — 2Ek), and 
their half-sum in the denominator. Since ujf = Ui — 2Q q , only one of the three terms can 
vanish at any time, and this diagram is clearly less singular than the one with the triple 
resonance. 

In fact, there are also higher order diagrams which have internal magnon lines or contain 
four-fermion vertices in low orders. They may formally be even more singular than the ones 
above. However, just as in the nonresonant section, we can omit them if we assume the 
large-S 1 limit, which was implicit in the considerations just outlined. 

B. Triple resonance 

In this subsection, we perform a detailed analysis of the diagram of Fig. ||(a), which gives 
rise to a triple resonance. The analytical expression for the diagram was given in Eq. fl5"T|). 

It may be useful to start with a remark on the related calculations of the two-phonon spec- 
tra in semiconductors. In that case the energy scale of the phonons is quite small compared 
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to the bandwidth in semiconductors. Therefore an expansion of electronic band structure to 
quadratic order near the minima is a priori warranted, and the resulting integrals are often 
doable analytically!!!!! 

In our case, however, the wavevectors of the magnons contributing to the two-magnon 
peak are of the same order as the Brillouin zone itself, and their energy is comparable to 
the fermionic bandwidth. Thus, we are forced to deal with a full spin density wave band 
structure without necessarily being able to expand it near some point. This complicates the 
integration considerably, although we will eventually find that some analytical results are 
still possible. 

We now turn to the general analysis of the possible singularities in Eq. (|51|). 
As we mentioned earlier, the three denominators in Eq. ([51]) vanish simultaneously when 
the following two conditions hold, 

Ui = 2E k ; Uf = Ui — 2Q q = 2E k+q . (55) 

Let us call k the value of k that solves the above equations, and expand the denominators 

to linear order about k . To study the conditions for a resonance, we set the numerator to 

a constant and obtain 

(4) f 1 
Mr ~ I £k {2w k0 ■ k - z5)((v feo + v fc0 _ g ) ■ k - i5)(2v fc0 _ 9 -k-i5) (56) 

where v ko = ^Hfc and v ko _ q = ^^\k ~q are the velocities at the two points in the momen- 
tum space where the denominators vanish. 

We now immediately see that this integral vanishes unless the two velocities are strictly 
antiparallel to each other. Indeed, let Vfc be in the x-direction. Then the first term in the 
denominator only depends on k x . If v ko _ q has a y-component, then only the second and 
third terms in the denominator will contain k y , and both with the same sign. In that case, 
if we do the k y integration first, the integral will vanish since the poles from both terms will 
lie in the same half-plane. 

The condition for the velocities at ko and (ko — q) to be antiparallel has a simple physical 
interpretation!! in the semiclassical picture, where both momenta and coordinates can be 
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assigned to the quasiparticle wavepackets. Suppose the incoming photon creates a particle- 
hole pair at k = ±k . The energies of particles and holes have different signs, so that the 
resulting electron and hole run apart with velocities +Vfc and — Vfc , respectively. After the 
emission of two magnons, with the momenta ±q, the electron and hole velocities become 
±Vfc _ g . Only if the electron and hole are "aimed back" at the point of origin, will they 
recombine to emit a final photon. Thus Vfc _ g has to be antiparallel to Vj. . 

Let us now count the variables in our equations for the resonance. For fixed ioi and 
ujf we have two resonance conditions, an additional condition that the two velocities be 
antiparallel, and also the energy conservation condition 2Q q = u = Ui — ujf. Thus, we have 
four equations for four unknown components of k and q. However, at a fixed u>i this system 
of four equations does not necessarily have a solution for an arbitrary ujf. For example, for 
the rotationally invariant magnon spectrum and parabolic fermionic bands, elementary con- 
siderations show that the two conditions in Eq. (|55|) plus the energy conservation condition 
fully determine |k |, |q|, and the angle between them, so for a fixed the condition on 
velocities has a solution only at a single Uf, and for that one Uf all of the magnons are at 
resonance at once. In the general case, the magnon spectrum is not rotationally invariant, 
but the triple resonance will still occur only in some range of the final photon energies Uf 
(and thus only for some magnon energies Q q ). We will find that this range is numerically 
rather narrow in most cases of interest. 

We now study the general solution for triple resonance. Consider first the mean- 
field (S = oo) forms of the fermion and magnon dispersions. We expand 2E/, = 



2y A 2 + At 2 (cosk x + cosk y ) 2 ~ 2A + 4JS(cosk x + cosk y ) 2 and introduce the reduced vari- 
ables 



Ui - 2A 
AJS 



oj f - 2A 
US 



(57) 



which count the photon energies from the band edge. Then Eq. (^) becomes 




(58) 
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where k = k , and where we introduced 1 = k — q and also restricted ourselves to the 
reduced Brillouin zone. 

The condition that the velocities ^Mr and ^ be antiparallel is 



sin k x sin l y — sin l x sin k y = 0, (59) 

provided that if k is in the first quadrant, then 1 is in the third (this makes the velocities 
antiparallel). 

Finally, the condition E\. — Ei = 2Vtk-i takes the form 



2 y 1 — ^(cos(fc x - l x ) + cos(k y - l y )) 2 = \ i -\ f . (60) 

These four coupled nonlinear equations is the set to deal with. There exists no general 
method for solving systems of nonlinear equations. However, in our particular case, much 
information can be obtained analytically by studying two high symmetry directions, which 
we do next. 

Let us first consider the case where q x — q y . It is easy to check that for such q, the solution 
for the resonance has a general form k x = k y (0 < k x < n/2) and l x = l y (— n/2 < l x < 0). 
The antiparallel velocities condition is automatically satisfied, and the remaining conditions 
yield 



cos k x — ^ ) cos lx — r} 1 
2^1 - cos 2 {k x - l x ) = 2sin(^ - l x ) = \ i -\ f . (61) 



These three equations trivially lead to 



Ai - 2 + y/3Xi(A - Aj) 



(62) 



where as before £\ = jfg is a reduced magnon frequency. 

By definition, the upper limit on \ is equal to 4. The lower limit is found by locating the 
point where A/ = A« — 2Q q becomes zero, since at that point one of the velocities vanishes, 
and after that we have parallel rather than antiparallel velocities and thus no contribution 
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from the triple resonance to the Raman vertex. Solving Aj = 2Q q together with Eq. fl62|) we 



find the lower boundary is at fl q — 1/2, Aj = 1. Eq. ( |5^) is graphically presented in Fig. [10. 

Another symmetry direction for which the condition on the velocities is satisfied auto- 
matically is that for which q is along the x or y-axis. In this case we can set k y = l y = q y = 0, 
and obtain a set of equations 

Aj = (1 + cos k x ) 2 , A/ = (1 + cos/ x ) 2 , 
Xi-Xf = 2^1-^(1 + cos(k x - l x )f. (63) 

It actually does lead to an analytic answer, but as it is rather cumbersome, we refrain from 
presenting it and consider instead only two special points. The first one is where the magnon 
energy is maximal, i.e., Aj — A/ = 2. This implies k x — l x = tt, and hence Aj = 9/4. The 
second special point is one at which A/ = 0, i.e., the termination point on the lower end. At 
that point l x = — n, and we easily obtain from Eq. fl63"l) that 

Aj(Aj + l) 2 = 16, (64) 

with the solution A, ~ 1.9. 



The whole curved tu g = /(A,) for q x = is plotted in Fig. [TO. 

Another part that is worth deriving analytically is the locus of the lower and upper 
termination points of the curves for different directions of q. It is easy to demonstrate that 
for the largest possible Aj = 4, the solution for any direction is Xf = 3, so that the curves for 
all directions of q meet at a single point when Aj = 4. To see this, note that Aj = 4 implies 
k x = k y = 0. After substituting this into Eqs. (|58D , (|59|), and (|60D, we obtain the claimed 
result. 

The lower termination point is Xf = 0, because the triple resonance clearly requires both 
the initial and final photon frequencies to lie above the gap. This condition corresponds 
to a straight line given by Aj = 2u = (oji — u)f)/4JS in the (Aj,a)) plane (recall that at 
resonance, cu = Cl q ). This by itself would allow < Aj < 2. A more detailed study shows, 
however, that the solution only exists for Aj > 1. Notice that curves for different magnon 
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directions terminate at different points on the line Aj = 2u>, and near Aj = 1 only a small 



fraction of magnon directions, namely q x ~ q y , allow a triple resonance there (see Fig. p~Q| ) . 
At the other end of this boundary (at Aj = 2), the numerator in fl5T|) vanishes. For all these 
reasons, Raman intensity is expected to be small in the region of low Aj, and it will not play 
any role in our subsequent analysis. 

Finally, we can also locate the region where our set of equations has a solution for the 
largest possible magnon frequency uj = Q q = 1. An inspection of Eqs. (|58|) , (fj9|), and (|60|) 
shows that the solution exists for 2 < Aj < 3. 

Combining all these analytical results, we obtain the allowed region of triple resonances 
in the (Aj, 2(D) plane. This region is shaded in Fig. [1^. Notice that although for 2 < 
Aj < 3 the solid line in the figure (the solution for q x = 0) practically coincides with the 
horizontal line 2lD = 2Q q = 2, this solid line is actually located below the maximum two- 
magnon frequency everywhere except Aj = 9/4. As an independent check, we also solved 
our equations numerically for a number of directions of q and found that the solutions were 
within the shaded region. 

The above solution for the triple resonance was obtained for the mean-field forms of 
fermionic and magnon dispersions. As we already discussed, quantum corrections will cer- 
tainly change the overall scales in the dispersions and thus modify the resonance conditions. 
However, we do not expect these modifications to be substantial especially near the upper 
termination point Aj = 4 as still, at maximum Aj we have k x = k y = 0, and hence curves 
with different directions of q will terminate at the same Aj. This region near maximum Aj 
will play central role in our subsequent analysis. 

Finally, it is worth noticing that the divergence we have found in Raman matrix element 
is an artifact of neglecting the damping of quasiparticles. If the damping was included, the 
divergence would be gone, and we would obtain instead only the enhancement of Mr in the 
shaded region in Fig. [10]. For this reason, later in the paper we will refer to an enhancement 
rather than to a singularity in Mr. It is important, however, that up to some damping levels 
the enhancement of Mr is going to be substantial, and the singularity analysis is a useful 
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guide to the actual physics. 



C. Two-magnon scattering 

So far we have considered the behavior of the Raman matrix element as a function of 
incoming frequency, neglecting finite state interaction, and selected the diagram which gives 
rise to a strong enhancement of Mr in some range of outgoing photon frequencies for a 



given u)i. At the same time, as we discussed in Sec. |III B| , the density of states is divergent 
at the magnetic zone boundary, and in the absence of finite state interaction this naturally 
gives rise to a divergence of R at to = 8JS. For S = 1/2, which we will assume below, this 
divergence is at u — 4J. 

If we now switch on the magnon-magnon interactions, the divergence at 4J will be 
gone, and we will obtain instead the maximum of R at a smaller frequency. In principle, 
this frequency should be different from that obtained in Sec. [1.1 FBI, as the Raman matrix 
element in Eq. ([51]) has a complicated dependence on the magnon momentum q, and the 
terms in the resulting series problem cannot be factorized. On the other hand, the location 
of the two-magnon peak at about 3 J is in agreement with the experiment aliS as well 
as numerical data.!9^ If so, the influence of the q-dependence of Mr on the magnon- 
magnon scattering should not be substantial. There will be also a feedback effect from the 
magnon-magnon scattering on the location of the region where the Raman matrix element 
is enhanced due to triple resonance. However, this feedback effect will be small in the 



semiclassical (large S) approximation, which, as we discussed in Sec. |III B| , is likely to work 
well even for S = 1/2. For all these reasons, for the rest of our discussion we adopt a 
semiphenomenological approach and assume that at any given Ui the Raman spectrum 
R(u) has two independent peaks: one is due to the triple resonance, which is particularly 
relevant to the first diagram in Fig. [| which has no finite state interaction, and the other, 
which for definiteness we assume to be at uj = 2.8J e //, is due to the magnon-magnon 
scattering (all other diagrams in Fig. ^). In essence, the approximation we made implies that 
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we are considering an "effective" Loudon- Fleury- like model, in which the Raman vertex is 
considered as constant in the calculations of the two-magnon peak position, but the overall 
magnitude of the vertex does indeed strongly depend on how close we are to the triple 
resonance region. 

Without considering in detail the effects of damping, we cannot conclude which of the two 
peaks is stronger. The experiments seem to indicate that the peak at 2.8 J e // is stronger, and 
the enhancement of the Raman matrix element can only be responsible for the asymmetric 
"shoulder-like" behavior of the two-magnon profile (see the discussion in Sec. [Vj). At the 
same time, if we fix u at 2.8J e //, as in Fig. [| and consider the variation of the two-magnon 
peak intensity as a function of the incident photon frequency Ui (or Aj), this intensity will 
clearly have a maximum where an uo = const line intersects the region of the triple resonance. 

To see where the maximum occurs, we neglect for simplicity the renormalization of J 
due to quantum fluctuations and draw a horizontal line in Fig. [10] at the reduced frequency 
uj = 0.7 (we checked that the results below are practically insensitive to whether we use J 
or J e ff in the magnon spectrum; however, if we used J e ff, we would also have to consider 
the 1/S renormalization for the fermionic spectrum, which does not lead to new physics but 
substantially complicates the calculations). We see from the figure that the line at ui = 0.7 
intersects the region of triple resonances in two places. One occurs at Aj very close to 4, i.e., 
when the incident photon is near the top of the fermionic band. The other occurs for Aj close 
to 1, i.e., not far from the bottom of the fermionic band. We know from the discussion given 
above that the region 1 < Aj < 2, especially for Aj close to 1, is rather "esoteric" in that 
only a small fraction of magnon directions allow for a triple resonance there (see Fig. |i~0| ). 



On the other hand, for A» close to 4, we will encounter no such problem, and we thus expect 
a much larger enhancement there. Thus, we focus entirely on the latter value of Aj ~ 4. 

We now consider how the two-magnon peak amplitude increases as Ui approaches the 
critical value where the triple resonance and the magnon-magnon peak occur at the same 
value of uo. First, since we are close to the top of the fermionic band, we can perform a 



quadratic expansion of the band structure near the maximum. Second, we see from Fig. |10 

33 



that near Aj = 4, the triple resonance region is very narrow, so that to first approximation 
we can consider that the triple resonance occurs only along a single line in the (Aj,u)) 
plane. As we discussed in Sec. |1V this would have been the case if the magnon spectrum 
was rotationally invariant, and one can easily check that it approximately is for Cl q ~ 0.7. 
Accordingly, we expect it to be a good approximation if we linearize the magnon spectrum 
around q = 0, 

n q ~ X. (65) 

The integral for the triple resonance can now be done analytically. The details of the 
calculations are presented in the Appendix. Here we quote the results. The resonance is 
allowed at the line 

A t = 4-(c-i) 2 , (66) 

which, in accordance with Fig. terminates at Aj = 4 for uj = 1/2, and intersects the two- 
magnon peak position (p = 0.7) at Aj = A 4 res = 3.96. The maximum of the two-magnon peak 
amplitude thus corresponds to the incoming photon frequency a;, that puts the particle-hole 
pair members very close to the tops of their respective bands. 

Let us now fix the two-magnon frequency at the two-magnon peak position and vary Aj 
(as in Fig. §).From the results in the Appendix, the Raman matrix element is then 

4 — A- 

M B ~ — Cl (67) 

R (Ai-Ar + ^) 3/2 1 ' 

The factor in the numerator comes from the vanishing of the term ^ in the numerator of 
Eq. (|51~D at k = 0, which is the value of k of the resonant particle-hole pair created by the 
incoming photon for Aj = 4, and also due to the fact that \ = 4 corresponds to fl q — \ = 
in which case the overall factor (1 — -^f ) = (1 — in Eq. ( |9~2"D in the Appendix vanishes. 

For Aj very close to A[ es , the term in the numerator can be considered as a constant, and 
we obtain 

M * ~ [A, -A,~ + «]»/» • (68) 
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In practice, however, the damping of fermions prevents the true singularity at Aj = A[ es , and 
the increase of M R can be measured at some distance away from A[ es . In this situation, the 
difference between X^ es and 4 can be neglected, and we obtain from ( |67D 

Mr m , (69) 

The sum total of those arguments is that the Raman intensity has at least an inverse linear 
singularity from triple resonance. 

R~ [Xi-Xr + t5]-\ (70) 

It is important to note that the usual single (Loudon-Fleury) resonance in 2D can at best 
give an inverse linear singularity of R for the spin density wave band structure, and only 



at the mean field level, due to the degeneracy of E k = yA 2 + e\. Once the mean-field 
degeneracy is broken by a better approximation, the single resonance is further reduced to a 
logarithmic singularity. Therefore the triple resonance singularity is expected to be at least 
as strong, and most likely stronger than any at the lower band edge. 

This is to be contrasted with the behavior of the optical conductivity ct(uj) (see Eq. ( |2"2"D 
above and Fig. |TT](a)). Unlike the Raman cross-section, it vanishes at the upper band edge 
Aj = 4. Most of the opitical weight in a(Aj) is concentrated towards the lower band edge 
Aj = 0, where it has a singularity in the absence of fermionic damping. Therefore, the optical 
conductivity resonates at the lower band edge while the Raman cross-section is strongest at 
the upper band edge of the coherent quasiparticle spectrum. 

Notice also that while the Raman intensity is symmetric with respect to the sign of 
A« — A^ es , the expression fl6"7|) for Mr, which includes only the leading singularity, is, strictly 
speaking, valid only in some region around A[ es . Away from this region, there are no reasons 
to expect the Raman intensity to be symmetric. On general grounds, we expect that the 
intensity above A[ es must fall down faster than for the same deviation from X\ es but at 
Aj < X^ es simply because no triple resonance is possible for Aj > 4. 

We now discuss in more detail how these (and other) results are related to experiment. 
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V. DISCUSSION 



We first summarize the main results of our work. 

We developed a general diagrammatic approach to Raman scattering in antiferromag- 
netic insulators. For photon frequencies small compared to the gap between conduction 
and valence band, we rederived diagrammatically the Loudon-Fleury Hamiltonian,0 which 
had been first derived from the Hubbard model by Shastry and Shraimanil in a different 
formalism. We also considered the location and shape of the two-magnon peak by studying 
the magnon-magnon scattering in a systematic 1/S expansion. The results of this study are 
consistent with the earlier considerations by Canali and Girvin.@ 

We studied for the first time the two-magnon Raman scattering in the so-called resonant 
regime when the incident and final photon frequencies are only O(J) apart from the Hub- 
bard gap. In this situation, the diagrams which are subleading in the nonresonant region 
become dominant. We identified the diagram which gave a dominant contribution to the 
Raman vertex in the resonant regime, and found the region in the (ui, uji — Uf) plane where 
the Raman vertex is divergent, unless one includes into consideration the damping of quasi- 
particles. The divergence is due to the simultaneous vanishing of all three denominators in 
the expression for the Raman matrix element. This phenomenon is called triple resonance. 
In the presence of a moderate quasiparticle damping, the divergence in the Raman vertex 
is gone, but the enhancement due to the triple resonance can indeed survive. 

Based on the experimental data, we then assumed that the location of the two-magnon 
peak in the resonant regime remains practically in the same place as at small photon fre- 
quencies, i.e., at Ui — ujf ~ 3 J. For nearly all incoming photon frequencies u>i falling within 
the Hubbard bands (i.e., corresponding to the features in the optical absorption spectrum), 
this magnon frequency is located below the region where our calculations show a triple 
resonance is allowed (see Fig. (0)). Based on these facts, we formulated an "effective" 
Loudon-Fleury-like theory in which the position and shape of the Raman profile in the two- 
magnon peak area of ~ 3000 cm (but not at higher magnon frequencies) is solely due to the 
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final state magnon-magnon interactions, while the triple resonance in the Raman vertex at 
larger u ~ 4000 cm -1 leads to the "shoulder-like" behavior in the magnon profile. When the 
incident photon frequency approaches u;™ aa: = 2A + 8 J, which corresponds to the particle 
and the hole being excited to nearly the tops of their respective bands, the triple resonance 
region allowed at a given Ui shifts down to lower transferred frequencies, and finally inter- 
sects the value of to where the two-magnon peak is located. We found that this happens 
when uJi/2 is very close to the top of the quasiparticle band, Ui = u^ es rs oj™ ax . Right at the 
crossing point, the intensity of the two-magnon peak amplitude diverges (in the absence of 
quasiparticle damping). 

As uJi approaches this critical value from below, the Raman intensity increases as R ~ 
(c<j[ es — Ui)^ 1 at some distance away from the resonance, and as R ~ {uj^ es — u>i)~ 3 in the 
immediate vicinity of the resonance. In a real situation, the true divergence will be gone 
due to damping of quasiparticles, and the two-magnon peak intensity as a function of uj{ 
will instead have a maximum at <jj\ es . We therefore do not expect to see cubic dependence, 
but inverse linear increase of R should be observable. 

We now analyze how these results are related to the experimental measurements in 
Fig. [I] and Fig. 0. In the Introduction we listed the key experimental features that required 
explanation. Here we list them again and comment on each of them. 

In Fig. |: 

(a) Asymmetry of the two-magnon peak profile: our theory predicts that for uo smaller than 
io\ e& the two-magnon peak profile should be asymmetric with a "shoulder-like" behavior 
at frequencies close to u — 4J, due to the triple resonance in the Raman vertex. This is 
consistent with experimental observations. In particular, the experimentally measured two- 
magnon profile in Pr 2 CuOi was analyzed^ and found to contain two peaks, a two-magnon 
peak at 3000 cm -1 , and a smaller one at 4000 cm -1 , which is precisely as expected from our 
calculations. Notice, however, that the calculation of the relative intensity of the peaks is 
beyond the scope of the present approach. 

(b) Selection rules: The leading diagram in the resonance regime contributes to scattering 
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in both Big and A\ g geometries. The signals in both geometries have been observed in 
the experiments. Recall that the Loudon-Fleury theory predicts scattering only in the B lg 
geometry. 

(c) Stability of the two-magnon peak profile: Away from the resonance, the stability of the 
two-magnon peak profile is consistent with our theory, as the enhancement of the Raman 
vertex for 2 < Aj < 3.5 occurs near uj = 4J. As uj{ approaches t^[ es , the situation be- 
comes more complex, as one cannot separately consider the peak due to magnon-magnon 
interaction and the enhancement of the Raman vertex due to the triple resonance. The 
experimental data seem to indicate that the two-magnon peak profile changes little as Ui 
sweeps through the resonance. This is not entirely consistent with our scenario, as the triple 
resonance region shifts to lower magnon energies u as the photon frequency Ui approaches 
io\ e& . However, the magnon profile itself in the resonant region can only be obtained by solv- 
ing coupled equations for the enhancements due to the triple resonance and the final state 
magnon-magnon interactions. We have not yet been able to accomplish that, and we leave 
it for future work. We regard this as the most important unresolved issue of all that remain. 
Notice also that the observed distortion of the Raman profile at higher photon frequencies 
is qualitatively consistent with our main conclusion that once the photon frequency passes 
the top of the coherent quasiparticle band, one must consider a qualitatively different set 
of electronic states. Our speculation at present is that at such frequences the dominant 
contribution to Raman scattering comes from an incoherent, diffusive hole motion, which 
can more easily couple to multi-magnon excitations. 
In Fig. |: 

(a) A single peak: Our theory predicts a single maximum in the two-magnon peak intensity 
measured as a function of the incident photon frequency. On the contrary, on the basis of 
general considerations in Raman scattering,^ we might have expected two peaks, one at LU{ = 
2 A, and the other at Uf — 2 A. These are the so-called incoming and outgoing resonances. 
In the more traditional area of phonon Raman scattering in semiconductors,0@'i! the two 
peaks are not always separately resolved. However, given the large magnon energy scale in 
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the cuprates, we might have expected them to be resolvable. An additional bit of evidence 
that triple resonance is a plausible explanation is that, experimentally, the two-magnon peak 
is seen in the interval of incoming photon frequencies about 4000 cm" 1 in width, which for 
J ~ 1000 cm~ l is of the same width as the interval 2 < A, < 4 in our theory, where the 
triple resonance is possible for all directions of the magnon momenta q. 
(b) Peak location: Our theory predicts the maximum of the two-magnon peak intensity, 
measured as a function of u>i, right near the upper edge of the quasiparticle fermionic band. 
At the same time, the coherent quasiparticle contribution to the optical conductivity o-{oj) 
was shown to have most of its weight located near the lower band edge and vanish at the 
upper band edge (as schematically illustrated in Fig. |ll|(b)). This is consistent with the 
experimental Fig. which shows that Raman scattering is strongest right at the upper end 
of those features in the optical data that can be interpreted as particle-hole excitations 
between the lower and upper Hubbard bands. The location of the Raman maximum versus 
the features in o~{uj) has been one of the greatest experimental puzzles in this area, and its 
solution in our approach serves as a partial verification of the SDW dispersion relation for the 
carriers which, despite much theoretical work, has not been well-established experimentally 
in these materials. 



In Fig. [12] we fit the published experimental dataES from Fig. |2| on the peak intensity in 
Y B 2 Cu^Oq at room temperature by our Eq. ([?(]). We see that a fit to the predicted inverse 
linear dependence is reasonably good. The inverse linear dependence starts from Ui ~ 2.4 eV 
and extends nearly up to the resonance frequency u^ es »s 3.1 eV. The effects of fermionic 
damping are probably relevant only in the immediate vicinity of the resonance. Note that 
the resonance value obtained from the fit is larger than the peak location at uj\ es ~ 2.8 eV 
in Fig |2|. In fact, the experimental data at high frequencies have been reexamined in recent 
studies^. It turns out that the resonance at room temperature occurs at a higher frequency 
than has been previously believed, and in particular the the intensity at uji ~ 2.9 eV is found 
to be larger than at u>i ~ 2.8 eV, unlike what one sees in Fig. |^. The new data suggest the 
resonance frequency value may exceed 3 eV, which would make it even more consistent with 
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the fit to the theoretical dependence. 

These recent experiments!! also measured the peak intensity in YB 2 CusOe_x at low 
temperatures (T ~ 5K). We fitted the new data by the same inverse linear dependence 
and found even better agreement with the experimental data for Ui between 2.5 eV and 
3 eV. This fit yields about the same resonance frequency u^ es ~ 3.1 eV as the fit at room 
temperatures. We actually do not expect that the resonance values at room temperatures 
and at 5K are the same, but we notice that the optical studies in La<iCuO± have showni^ 
that the optical absorption peak shifts to higher frequencies by less than 0.1 eV between 
room temperatures and 122 K, and one can expect the peak in Raman intensity to shift 
by about the same amount. This small shift is indeed within the combined accuracy of our 
theory and the data available for the fits. 

Notice that our theory also predicts that the Raman intensity should have an additional, 
smaller, peak at uo ~ 2 A ~ 1.7 eV due to the Loudon-Fleury mechanism. There is no 
evidence for the peak at u>i = 2A in the room temperature data on Fig Q Preliminary 
results of recent low-temperature experiments^ indicate a possibility that there is in fact a 
second, smaller peak at u>i ~ 1.6 eV which is consistent with our prediction. 

There are, however, several experimental results which are beyond the scope of our 
approach. First is the width of the two-magnon peak, which exceeds the prediction of the 
spin-wave theory even if we take into account the resonance in the Raman vertex. As we 
already remarked above, simply invoking quantum fluctuations for spin- 1/2 is unlikely to 
improve the spin-wave theory, since direct numerical calculations of the Raman spectrum 
on finite clusters also predict a narrow peak in this case.0 Another experimental fact is the 
existence of a considerable Raman signal R(u) above the maximum possible two-magnon 
energy (i.e., above 4J e ff). Canali and GirvinS performed a very detailed study of the 
effects of four-magnon scattering within the context of the Heisenberg model for spins and 
the Loudon-Fleury coupling of spins to light. They found that this scattering can give rise 
to the Raman cross-section above 4J e //, but the intensity of the Raman signal was found 
to be too small to fully account for the experimental data. Thus one apparently has to go 
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beyond the minimal model. 

In the approach put forward by Singh,!!! further neighbor terms in the effective Hamil- 
tonian for coupling of light to spins are postulated, and the coefficients adjusted so that the 
moments of the Raman spectrum evaluated by the series expansion agree with observations. 
Such next-nearest-neighbor terms broaden the two-magnon peak and also automatically 
break the selection rule for the Loudon-Fleury Hamiltonian which disallows scattering in 
the Ai g configuration. Of course, a question then arises of the origin of such terms and 
their relative strength. Notice that the further-neighbor terms are in fact also effectively 
present in our momentum-space description, as it clearly follows from the Shraiman-Shastry 
formulation of the problem in which the non-Loudon-Fleury terms in the Raman vertex 
correspond to the effective spin interactions between further neighbors. 

In the bulk of the paper we "factorized" the problem by setting the Raman vertex to 
its Loudon-Fleury form in the calculations of two-magnon profile right near the peak. This, 
however, is only an approximation, and in view of Singh's results, we can expect some 
broadening of the two-magnon peak in the more sophisticated calculations along our lines, 
which would treat the resonance and the final state interactions together. Simultaneously, 
we might also expect a shift of the two-magnon peak position from 2.8J e ff. In the bulk of 
the paper we used 2.8J e // because of the agreement with neutron experiments. We note 
however, that the only observations to date of the two-magnon peak in cuprates were made 
at or near the resonant frequencies, which are in the visible light region in cuprates. In that 
sense, we may not even know what the "true", i.e., nonresonant peak shape and location 
are. 

The scattering in B\ g and A\ g geometries in the SDW technique has been studied by 
Kampf and Brenig.il They, however, focused on high transferred frequencies, comparable to 
the Mott-Hubbard gap, and did not include collective spin fluctuation modes (i.e., magnons) 
in their theory. In view of this, a comparison of our results with theirs is not possible. 

In a rather different spirit, it was observed in a paper by Weber and FordS that even a 
rather small magnon damping introduced phenomenologically into the equations broaden's 
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the Parkinson peak considerably. More recently, MerlinEa argued that phonons may be 
the source of that damping. Should the electron-phonon interaction prove relevant to the 
problem, this will only lend credence of the key phenomenological assumption of this work, 
namely that on a "first pass" at the resonant scattering problem, it is best not to try to 
deal simultaneously with two quite possibly distinct issues of the two-magnon peak shape 
and its strength variation with the incoming photon frequency. 

Finally, an experimental fact which needs answering is the existence of the Raman sig- 
nal at very high incident and transferred frequencies, of the order of 3.5 eV and 1 eV, 
respectively. The largest intensity at such high frequencies was obtained in the "chiral" 
A 2g channel. Khveshchenko and Wiegmann recently performed an elegant study of the 
contribution to Raman vertex in this frequency range from chiral spin fluctuations in the 
magnetically ordered phase.il For our considerations, it is essential that at such high fre- 
quencies, the incoming photon is above the coherent quasiparticle band (where our theory 
applies) and the motion of the carriers that subsequently couple to the magnons is itself 
incoherent in character. This and other interesting issues require further experimental data 
as well as an improved theoretical understanding. 



VI. ACKNOWLEDGEMENTS 

It is our pleasure to thank G. Blumberg, C. Canali, S. L. Cooper, D. Khveshchenko, 
M. V. Klein, R. Liu, R. Martin, R. Merlin, H. Monien, D. Pines, S. Sachdev, C. M. Varma, 
P. Wiegmann, and A. Zawadowski for useful discussions and comments. We are grateful 
to G. Blumberg, M. V. Klein, S. L. Cooper and R. Liu for providing experimental data 
for our Figs. |], and [TJ. D.F. is supported by the Texas Center for Superconductivity 



at the University of Houston. Part of the work has been done while A.C. was at Yale 
University, where he was supported by the National Science Foundation (NSF) Grant Nos. 
DMR-8857228 and DMR-9224290, and while D. F. was at the University of Illinois, where 
he was supported by the NSF (DMR 91-20000) through the Science and Technology Center 



42 



for Superconductivity at the University of Illinois at Urbana-Champaign. 



VII. APPENDIX 



In this Appendix, we obtain Eq. (p7|) for the Raman matrix element near the intersection 
between to = 2.8 J e ff, where the two-magnon peak occurs, and a (narrow) region in the (A i; to) 
plane where triple resonance is allowed. For simplicity, we will neglect the difference between 
J and J e ff- We see from Fig. |TD| that this intersection occurs very close to the top of the 
quasiparticle fermionic band E k = (A 2 + 4t 2 (cos k x + cos ky) 2 ) 1 ^ 2 , i.e., near k = 0. Once this 
is established for this band structure, we can get a good approximation for the integrals by 
expanding the energies to quadratic order around k = 0. 

The integral of interest is Eq. (|5lD for the diagram in Fig. ||(a): 

M R = SiT' (t ' ^ ( ^ ' ^ " (71) 

^ (u)i - 2E k + iS) (u f - 2E k _ q + iS) (ui - Vt q - E k - E k _ q + iS) 

We have upon expanding 

e k = -2t(cos k x + cos k y ) ~ -At M - M , (72) 



and (we consider S = 1/2) 



k 



2 N 



E k = ^A 2 + e 2 ^A + U\ [ l- Y j. (73) 

We will expand both E k and E k+q , which implies that we also assume that the typical magnon 
momenta are not too large. This, as we discuss in the bulk of the paper, is consistent with our 
explicit result that the width of the region where triple resonance is allowed is very narrow, 
nearly a single line. To the same accuracy we can also linearize the magnon spectrum, 
f2q = \f2Jq. However, we found it convenient to keep using the general form Q q for the 
magnon frequency in some of the formulae. 

To study the singular behavior, we first set the numerator to a constant. Expanding in 
powers of k in ( [73| ) and introducing, as before, the reduced variables, Aj = (a;, — 2A)/2J, 
Xf = (ojf — 2A/)2J, and Cl q = Q q /2J, we obtain the following integral, 
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d 2 k- 



(Xi - 4 + 2k 2 + i8) (A/ - 4 + 2(k - q) 2 + i5) (A* - Q q - 4 + k 2 + (k - q) 2 + iS) 

(74) 

Somewhat similar integrals appear in the three-dimensional problem of two-phonon reso- 
nant Raman scattering, where the incoming photon frequency is tuned to match the critical 
points in the semiconductor band structure. The case when the hole band has an infinite 
mass was done analytically by R. Martin.il The case of equal effective masses in two di- 
mensions has not, to our knowledge, been studied analytically before, and we present the 
derivation in some detail. 

We first observe that one of the terms in the denominator of the integrand is half the 
sum of two others. The integrand can thus be put in the form 



1 



1 



1 1 



We now recall Feynman identity familiar from Quantum Electrodynamics, 

l r 1 . 2x 



(75) 



dx 

P 2 Q Jo (xP + (l-x)Qf 
Using this identity and completing the squares containing k, as necessary, we obtain 



d 2 k / dx 
Jo 
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L (A* - 4 + i5 - xn q - (f - x)q 2 + 2(k - f q) 2 

x 



'X, - 4 + iS - (2 - x)n q - (f - x)q> + 2(k - (1 - f )q) 2 
Shifting the variables of integration and using the integral 
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we obtain 
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Completing further the squares with respect to the variable of integration x, we rewrite 
the r.h.s. of fl75| ) as I — I\ + h, where 



7T 



7i = — / <ix 
4 jo 



.2; 



Ai-4+f +i<J-£(x-(l 



1 2 ' 



7T 



h — t / ^ x 
4 jo 



Ai-4+f (l-£) + 
At this stage, we have reduced I to two integrals of the general type 



J 



dx- 



x 



/o ""[C - F 2 {x- D) 2 } 2 ' 
In both 1 1 and I2, we have the same 



and F 2 = y- However, the TJ's are not the same. We have 



and 



Do. = 1 



in ii and ^2, respectively. 
We now rewrite J as 



J 



dx- 



x-D d r 1 D 

^~dCJo C-F 2 {x-D) 2 ' 
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Performing the integrations and partial differentiation with respect to parameter, we obtain 
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2F 2 (C - F 2 (x - D) 2 ) 



1 D 
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x-D 



2C (C — F 2 (x — D) 2 ) 
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4Fcf \F(x-D) - VC, 







(86) 



The integral I is a sum of two such integrals and is thus 
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h + h 
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It turns out that the first four terms inside the braces add up to zero. To see this, we 
observe that the definitions of D\ and D 2 given above imply that 
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which in turn implies that (1 — Di) 2 = {I — D 2 ) 2 , and 

Di D 2 



+ 



We also use the identity 



2(1 -A) 2(1 -D 2 ) 



1 



+ 1 = 0. 
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whenever terms of this type occur. Using all these identities, we can easily show by a direct 
calculation the cancellation claimed above. 
We now have a compact result 
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Substituting the full expressions for , D 1 , D 2 , and F, we obtain 
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The vanishing of C corresponds to the triple resonance. Near this point, I behaves as 



46 



The expression above is incomplete until we specify the branches of the logarithms and 
the square roots. A careful examination of the derivation shows that we must choose the 
square root branch which is positive for C > and the branch of the logarithms which tends 
to zero when its argument goes to unity. The answer for other real values of C is obtained 
by setting C — > C + iS and then analytically continuing C along the real axis. 

We now turn to the form-factor in the numerator of Eq. ([711) . Clearly, for any X r { es < 4, 
this form-factor can be considered constant in the immediate vicinity of A[ es . At the same 
time, since the form-factor vanishes at the top of the quasiparticle band, it will contribute 
a small factor to Raman intensity for A^ es ~ 4. This will effectively change the singularity 
of the Raman matrix element at some distance away from the triple resonance, since for Aj 
away from 4 the difference between A[ es being exactly 4 or merely close to it is not noticeable. 

Near A[ es = 4, we set (fi q €k~ q — X q tk) = const, linearize ~ k and ^ q ~ k — q, 
expand the numerator terms to quadratic order, and obtain 

[k-e,][(k-q)-e>] 



M. 



R 



E 



t (Ai - 4 + 2k 2 + iS) (X f - 4 + 2(k - q) 2 + id) (\i - Q q - 4 + k 2 + (k - q) 2 + i5 

(93) 



in place of Eq. (ff^) above. Then we proceed using Feynman parametric representation as 
in the case with a constant numerator above. The two shifted k variables in Eq. ([FTP were 
k' = k — |q and k" = k — ( 1 — |J q. We can alternately express the numerator in the 
integral above in terms of these as 



[k-e^k-q)-^ 



/J 



and 



k' + fq 
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l-f]q 



3 / 
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1 - | ] q ] • e, 
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(95) 



When we now integrate over k' and k" in the B\ g scattering geometry, the angular 
integrals will only leave the terms — | (1 — | )(q • e^)(q • e^) in the numerator. Thus, we 
obtain the same integrals as in Eq. fl80D above, except that instead of x in the numerators 



47 



.2 



we will have — ^-(1 — |) . The angular factor involving magnon momenta, (q- e^) (q- e^) , when 
computed in the B\ g geometry, is just an expansion of the Loudon-Fleury vertex Mr 1 of 
Eq. (|39|) to quadratic order in q, which is consistent with the linearization of the form-factor 
done in the beginning of this calculation. 

The integrals over the parametric variable x can now be done by elemetary means similar 
to those employed above for a constant numerator case. The resulting expressions are, 
however, too lengthy to be worth presenting here. At the same time, as we discussed above, 
what is of real interest to us is the modification of the singular behavior of Mr at the values 
of A, at which one can neglect the difference between \ — A[ es and \ — 4. A simplest way to 
proceed is to consider the case A[ es = 4, when we also have Q q = \. In that case Di = in 
Eq. ([S3l), and D 2 = 2 in Eq. (|8~4f) . Thus, of the two integrals like the one of Eq. (^l|) above, 
but now with — \{1 — f ) instead of x in the numerator, only the one with D\ is singular 
for small C . The singularity is, furthermore, picked up near x = 0, so that we can also omit 
the factor (1 — |) in the numerator. We end up with 



with C = Xi - 4 + id. 

Performing an elemetary integration, and omitting an irrelevant nonsingular term, we 
obtain 



which has a leading square root singularity This is the result needed in the main text. 
Combining now Eqs. (0) and (p7|) and restricting to only a leading singularity, we obtain 
the result quoted in Eq. (|67D in the main text. 
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FIGURES 

FIG. 1. A typical Raman cross-section as a function of transferred photon frequency. A 



two-magnon peak is clearly seen. Data courtesy of the authors of Ref. 18. 



FIG. 2. The strength of the two-magnon peak as a function of incoming photon frequency. Also 
shown is the imaginary part of the dielectric constant. Data courtesy of the authors of Ref. |l8|. The 
data in the Figure were obtained at room temperatures. At low T, both the peak in the optical 
absorption and the Raman profile shift to higher frequencies by about O.leylElH The position 



of the peak in YBCO has been reexamined in recent studies 65, which place the maximum at 
uji ~ 3 eV rather than 2.8eU, as in the figure. 

FIG. 3. The real-space picture of the two-magnon scattering. Dashed lines denote photons, 
and wavy lines denote spin waves. The explanations are given in the text. This picture is valid 
away from resonance. 

FIG. 4. The diagrams for two-magnon emission which contribute to the Loudon-Fleury Hamil- 
tonian at small incident frequencies. Each fermion can belong to either the valence (dashed line) 
or conduction (solid line) band. The emitted magnons are denoted by the solid wavy lines, and 
the incoming (ojj) and outgoing (u>f) photons by the dash-dotted lines at the ends of the diagrams. 
Additional graphs are obtained from each diagram by the fermion flow reversal and/or flipping all 
the spin labels. 

FIG. 5. Additional graphs which are of the same order of magnitude in t/U as those of Fig. 3, 
but have extra smallness in 1/5. The four-fermion vertex is the Hubbard interaction term U. 
Graph (a) is obtained from Fig. 3(c) by the insertion of a particle- hole bubble in place of a single 
four-fermion interaction. A sequence of such bubbles sums up to an internal magnon line; an 
example is given in (b). 

FIG. 6. A series of diagrams with the final-state magnon-magnon interactions. The black dot 
is the Loudon-Fleury vertex for the photon-magnon interaction. 
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FIG. 7. The profile of the scattering cross-section in the B\ g geometry without the final-state 
magnon-magnon interactions (dashed line) and with final-state interaction (solid and dotted line) 
The solid and dashed lines were obtained, correspondingly, by using the full expression for I(u>) in 
(refl) and its expansion near the top of magnon band which is in line with 1/S expansion. Notice 
that while the position of the two-magnon peak is the same in both cases, the expansion near the 
top of magnon band yields much broader two-magnon peak 

FIG. 8. Diagrams which become important at resonance. All three denominators in the diagram 
(a) can vanish simultaneously which is known as a triple resonance. This diagram is dominant in 
the resonance region. Additional relevant diagrams are generated by the spin label reversal (in 
both (a) and (b)) and/or the emission of two magnons from the valence band fermion line (in the 
diagram (b)). 

FIG. 9. The lowest order vertex (a) and self-energy (b) corrections to the mean-field theory. 

FIG. 10. The triple resonance region (shaded) in the (Aj, (u>i — ujf)/4:JS) plane where 
Aj = {uji — 2A)/4J5. Dashed line - solution for a triple resonance for magnon momentum q x = q y ; 
solid line - solution for q x = or q y = 0. The horisontal line corresponds to the position of the 
two-magnon peak at {uji — u>f)/AJS = 2lu = 1.4. 

FIG. 11. (a) The plot of Eq. ( p2|) for the optical conductivity, in the limit t <^iU . It has a square 
root divergence at the lower band edge, and vanishes at the upper band edge. The divergence will 
be washed out in practice, but we can still expect most prominent features to be located towards 
the bottom of the SDW quasiparticle band, (b) A schematic illustration of the relative location 
of the strongest Raman and optical features as a function of the incident photon frequency. Their 
relative position is a key result of this work. The two solid lines represent the fermion dispersion 
curves in the conduction and valence bands. 
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FIG. 12. A fit of the experimental dependence of the inverse two-magnon peak intensity in 
YB2C113O6 from Fig. [2] to the theoretical l/(u;[ es — ui) dependence, Eq. (ffO|) . The value of w[ es 
from the fit is ~ 3.1eV. 
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